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Abstract 

We investigated theoretically the effect of covalent edge functionalization, with organic 
functional groups, on the electronic properties of graphene nanostructures and nano-junctions. 
Our analysis shows that functionalization can be designed to tune electron affinities and ion- 
ization potentials of graphene flakes, and to control the energy alignment of frontier orbitals 
in nanometer-wide graphene junctions. The stability of the proposed mechanism is discussed 
with respect to the functional groups, their number as well as the width of graphene nanos- 
tructures. The results of our work indicate that different level alignments can be obtained and 
engineered in order to realize stable all-graphene nanode vices. 
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The successful fabrication of the two-dimensional system graphene l and its excellent mechani- 
cal and electronic properties have aroused great interest in several different fields, both from funda- 
mental research and for technological applications.^! Quasi-one-dimensional graphene nanorib- 
bons (GNRs) with semiconducting behavior can also be obtained, 6 7 with electrical band gap tun- 
able within a wide range,^ and their optical properties, even if less intensively studied, offer 
stimulating perspectives. Q^El One main challenge is to gain control of the graphene edge mor- 
phology at the nanoscale, in view of designing all- graphene nano-devices. In this sense many 
recent improvements in production techniques were reported, including both top-down^ED an( j 
bottom-up^2E3 approaches, which point to the possibility of exploiting the exceptional properties 
of graphene nanostructures in actual devices (see e.g. RefsJ23HZ5}). 

In view of electronic and optoelectronic applications, it is important to identify further ap- 
proaches for an efficient engineering of the electronic states of a graphene nanostructure. It 
has recently been shown that the electronic properties of graphene-based systems can be con- 
veniently modulated by means of organic functionalization.22H22] Until now, a number of proce- 
dures have been proposed exploiting either surface functionalization through chemisorption of 
molecules,'20E2] which is expected however to damage the transport properties of the system, or 
non-covalent interactionJ^HS] which usually induces doping effects but does not lead to stable 
structures. A more effective strategy could be edge functionalization,'^^ which may offer promis- 
ing perspectives also in so far unexplored application fields, such as the production of GNR-based 
nanoj unctions. 

In semiconductor physics, heterojunctions are usually classified on the basis of the relative 
alignment of valence and conduction band states at the interface of the two materials. Heterostruc- 
tures where the band gap of one semiconductor is fully contained in the gap of the other are 
commonly indicated as type I or straddling. On the other hand, the so called type II or staggered 
alignment is achieved when the top of valence band and the bottom of conduction band of one 
semiconductor are both higher in energy compared to the corresponding states of the other, so that 
the top of the valence band of the junction lies on one side of the interface, and the bottom of 



conduction on the other side. In this case, electrons and holes on the frontier orbitals are spatially 
separated, being located on opposite sides of the junction, which may promote relevant processes, 
such as e.g. photoinduced charge separation, as requested by light harvesting applications, or 
photocatalysis. These concepts can be extended to graphene nanoj unctions, where the seamless 
covalent transition from two graphene nanostructures with different electronic characteristics can 
create, as in conventional semiconductor heterojunctions, different kind of alignments. Indeed, 
interfacing semiconducting GNRs with different widths, and consequently different energy gaps, 
provides a type I level line-up Q2EQI j n ^ s paper we propose a proof of concept for a viable method 
to realize all-graphene type II nanoj unctions. For this purpose we investigate the effects of stable 
edge functionalization through organic electron-donating and electron- withdrawing groups, which 
modify electron affinities and ionization potentials in comparison with the H-terminated GNR. The 
proposed model for chemical engineering carries a high potential for producing stable graphene 
nanostructures, preserving their flexibility and their intrinsic electronic features, sp 2 carbon net- 
work being unaltered upon functionalization. 

Methods 

As prototypes for graphene nanoribbons we consider sub-nanometer-wide graphene nanoflakes 
(GNFs), that is, finite systems, and we focus in particular on H-terminated GNFs Q a )] 5 organically 
functionalized GNFs (TJb)], as well as on the nanoj unction obtained by functionalizing one half of 
the flake with -COCH 3 groups [TRc)]. In this case, we will be interested in the frontier molecular 
orbitals (that stand here for the conduction and valence band states) and their localization and 
energy alignment across the differently terminated regions. In order to have access to the electronic 
energies, we performed calculations in the frame of the semi-empirical AMI schemeJ 41 ^ AMI 
is a well tested method based on the Hartree-Fock scheme, 43 suitable for finite structures ;S3H25] 
it allows one to evaluate the ground state properties of the system and to compute the electrical 
band gap Eq as the difference between electron affinity (EA) and ionization potential (IP). The 



latter quantities are expressed in terms of the total energy of the neutral [E(0)] and charged systems 
[E(±l)], in the usual way EA = E(0) - E(-l) and IP = E(+l) - E(0). For all the considered structures, 
we performed full geometrical optimization to evaluate E(0), with a threshold for the forces of 0.4 
kcal- mol^VA. The E(±l) are obtained at the fixed atomic positions optimized in the ground state, 
i.e. vertical EA and IP are computed. 



12 3 4 5 




\ 


t: 



Figure 1: (Color online) (a) Hydrogenated C 80 H 28 graphene nano flake (GNF); the numbers iden- 
tify the hydrogen atoms that can be substituted by the organic groups upon functionalization. 
(b) Fully functionalized GNF obtained by replacing each second H atom terminating the edge 
[i.e. 1-5 and l'-5' in panel (a)] with a -COCH 3 radical; in the inset the atoms of the radical are 
highlighted, (c) Nanoj unction obtained by fully functionalizing one side of a C 12 2H 40 GNF with 
-COCH 3 groups. 



Results and Discussion 



Our reference system is the C 80 H 28 flake l[TVa)], characterized by armchair- shaped H-saturated 
edges along the length (x direction in [T]). We checked that the chosen length/width ratio is such 
that the zigzag-shaped end-borders (y direction) do not affect the intrinsic properties related to 
the armchair edges along the length. Thus, C 80 H 28 is a suitable model for an armchair GNR of 
width parameter N=7, N being the number of dimer lines across the ribbon width, according to the 
standard notation.^ It corresponds to an effective width of ~ 7 A, defined as the distance between 
C atoms on opposite edges, belonging to the same zigzag chain along y. Next, we addressed 
the same GNF upon covalent edge functionalization with organic groups. We started focusing 



on molecules containing a ketone functional group, characterized by high electronegativity of its 

carbonyl oxygen atoms, and in particular we considered the methyl ketone group -COCH 3 [inset in 

[jjb)], which is known to be more stable than the simpler formaldehyde group -COH. We simulated 

full functionalization of the C 80 H 28 by replacing each second H atom along the length with methyl 

ketone, corresponding in this case to ten anchored groups per flake [C 80 H 18 (COCH 3 ) 10 ], i.e. five 

per edge side, as shown in[]Jb). We computed Eq, EA and IP, and we compared these quantities 

to those evaluated for C 80 H 28 by calculating AEA (AIP) as the difference between EA (IP) of the 

functionalized and hydrogenated GNF, respectively. A positive shift of both EA and IP is revealed, 

as shown in the second column of \\\ 

Table 1: (Color online) Energy gap (Eg) and variations in EA and IP (AEA and AIP), 
computed with respect to the reference hydrogenated C 80 H 28 (H), of GNFs fully function- 
alized with different organic groups. Among the considered groups, -COCH 3 , -CH=CF 2 , 
-OCF=CF 2 and -OCH 3 give rise to a type II offset, with a positive (negative) sign of AEA 
and AIP for electron-withdrawing (-donating) groups; -CH=CH 2 and -(CH) 3 =CH 2 produce 
a type I offset. 
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In addition to the methyl ketone group, we considered the case of other organic functionaliza- 
tions. We focused on few prototypical substituents, such as (i) a methoxy group (-OCH 3 ), charac- 
terized by an O-bridge ester bond, which is known to be electron-donating with respect to aromatic 
molecules ;EQ(ii) a -CH=CF 2 group, which is expected to be electron- withdrawing, due to the high 
electronegativity of its F-termination; (iii) a -OCF=CF 2 group, which combines the presence of 
O-bridge ester bond and halogen termination. Due to their characteristics, we expect group (i) to 
provide a similar effect compared to methyl ketone, i.e. to rise EA and IP with respect to the hydro- 
genated flake, and group (ii) to do the opposite. In the case of -OCF=CF 2 the coexistence of both 



O-bridge bond and F-termination makes the picture less predictable. We also consider -CH=CH 2 
and -(CH) 3 =CH 2 functionalizations which do not introduce any foreign species in the system and 
preserve its C-sp 2 conjugation through the alternation of single and double C-C bonds. For each 
chosen group we addressed full functionalization, analogously to the case of -COCH 3 illustrated 
in[jjb), and we computed Eq, EA and IP. 

By inspecting the complete results in[T] two main features arise. The first is related to a decrease 
of the energy gap, independently of the anchored group: this is ascribed to an increase of the 
effective width of the flake, compared to the fully hydrogenated case. The other important outcome 
concerns the sign of AEA and AIP: we observe positive (negative) values of both AEA and AIP, 
corresponding to a downshift (upshift) of both HOMO and LUMO for groups with a predominant 
electron-withdrawing (-donating) character, as well as for -OCF=CF 2 . On the contrary, -CH=CH 2 
and -(CH) 3 =CH 2 functionalizations produce AEA and AIP having opposite signs: given the simple 
bond- alternating character of the functional group, the major result is related to the extension of 
the wave function to a larger effective width, which thus reduces Eq in agreement with results for 
one-dimensional GNRs.^l 

The analysis of Mulliken charge population shows that a charge redistribution at the edge oc- 
curs upon functionalization. In particular Mulliken charges of functionalized C edge atoms become 
more positive, being these atoms bonded to more electronegative species than H. Larger differences 
(> 0.2 e~) are observed when functional groups are bonded to the flake through an O-bridge ester 
bond, i.e. -OCH 3 and -OCF=CF 2 . On the other hand, moieties bonded through C atoms induce 
slighter Mulliken charge redistribution, of the order or less than 0.1 e~ . It is important to highlight 
that these redistribution effects vanish within one aromatic ring. 

We also performed some tests on functionalized C 80 H 28 at decreasing number of functional 
groups (see [2]). We focus here in -COCH 3 as an illustrative example but this behavior is indepen- 
dent of the anchored group. We take the hydrogenated flake and, retaining the symmetry of the 
system, we introduce increasing even numbers of functional groups, up to ten, which corresponds 
to full functionalization for the chosen C 80 H 28 flake (seeHTb)). We start from two anchored groups 



Table 2: Energy gap (E G ) and variations in EA and IP (AEA and AIP) of the GNF function- 
alized with an increasing number of -COCH 3 groups, up to 10 (full functionalization), with 
respect to hydrogenated C 80 H 28 (H). 
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5.13 
0.00 
0.00 


5.01 
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in positions 3 and 3', after the scheme inPiTa), which corresponds to a -COCH 3 substituent each 
fifth H atom. We then introduce four and six groups in positions (2, 2', 4, 4') and (1, 1', 3, 3', 
5, 5') respectively. In the case of eight bonded groups, in order to maintain the symmetry of the 
graphene flake, each H atom in positions (1, 1', 2, 2', 4, 4', 5, 5) is replaced by a -COCH 3 . The 
results, reported in[2j show a monotonic growth of both AEA and AIP with increasing number of 
covalently bonded groups, with an almost linear dependence for AIP. At the same time, we no- 
tice the monotonic decrease of Eq, consistently with the above mentioned picture related to the 
increased effective width of the GNF. 

We also discuss the dependence of the electronic properties of the GNF, specifically related to 
EA and IP, by directly increasing its width. For this purpose, we performed a systematic investi- 
gation on GNFs, characterized by fixed length (along x direction in[T|) and increasing width (along 
y direction), fully functionalized with 16 -COCH 3 groups.— We studied eight GNFs ranging from 
width parameter N=4 (less than 4 A wide) to N=l 1 (about 12 A wide) and we computed EA and IP 
for each of them. As shown in [2j the calculated values for EA and IP have an oscillating behavior 
characterized by the modulo 3 periodicity already encountered in infinite GNRs EE3|54] ^^ ^g 
smallest value of the energy gap pertaining to the N=3p + 2 family (p integer). An approximate 
~ 1/vv dependence (being w the GNF width) is recognizable for the energy gap, consistently with 
well known theoretical^ and experimental^ results. The values of AEA and AIP, represented in 
[2]by the spacing between the respective marks, show a faster decay, which may be compatible with 
a local dipole mechanism driving the offset between functionalized and hydrogenated flakes. For 
instance, for the 3p + 2 family, we get values of AEA=0.85 eV and AIP=0.67 eV for the smallest 
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Figure 2: (Color online) Negative of electron affinity and ionization potential (-EA and -IP, re- 
spectively) of hydrogenated (H) and (-COCH 3 ) functionalized GNFs of increasing width. We also 
indicate the parameter N, i.e. the number of C-C dimers along the width. 



computed flake (N=5, w=4.9 A) and AEA=0.53 eV and AIP=0.43 eV for the largest one (N=ll, 
w=12.3 A). From the results in|2| we can extrapolate non negligible effects of covalent function- 
alization also for nanometer- sized graphene flakes, which can be considered realistic for current 
nanofabrication approaches. 

Finally, we directly investigate the type II nanojunction obtained by connecting a fully hydro- 
genated GNF to a fully functionalized one, having the same width (see[TJc)). For this purpose we 
adopted a flake of the same width of C 80 H 28 but increased global length, in order to ensure that each 
side of the junction is long enough to confidently reproduce the bulk properties of the flake (i.e. 
same single-particle orbitals in the energy range close to the energy gap). The resulting system ac- 
commodates full functionalization of one half flake with eight radicals, four per side (see[TJc)). We 
first performed full geometrical optimization and then focused on the ground state electronic prop- 
erties of the junction. In particular we are interested in the energy and wave function localization 
of the frontier states, as summarized in[3j The analysis of the wave function distributions clearly 
reflects the type II character of this junction. In fact we observe that the top valence states, HOMO 
and HOMO-1, are localized on the hydrogenated side of the junction, while the bottom conduction 
states, LUMO and LUMO+1, are localized on the functionalized side. On the other hand the next 
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Figure 3: (Color online) Energy levels (center) and wave function distributions (sides) of the fron- 
tier states of the nanojunction (]Jc)]. The dotted grey line in the energy level diagram indicates 
the reference (from the homogeneous GNFs) for HOMO and LUMO of the hydrogenated (H) and 
functionalized (-COCH 3 ) sides of the GNF. The energy positions of the six frontier orbitals repre- 
sented at the side panels are marked by the heavier (blue) lines: the solid intervals indicate finite 
wave function probability, while dashed intervals indicate vanishing wave function probability. 

states HOMO-2 and LUMO+2 show already resonant character. This feature can be explained 
in terms of the absolute energy levels of the bulk of the two sides, as indicated in the graph of 
[3] we represent by the dotted grey line the gap of the homogeneous systems, hydrogenated and 
functionalized, evaluated from the HOMO and LUMO levels of the fully hydrogenated [C 116 H 40 ] 
and functionalized [C 116 H 24 (COCH 3 ) 16 ] GNFs. On top of this we plot the level spectrum for the 
junction-GNF, highlighting the localization character of the related state through blue solid (finite 
probability) or dashed (zero probability) lines. We see that HOMO and HOMO-1 lie in the for- 
bidden gap for the functionalized side and are consistently localized on the allowed hydrogenated 
half; the opposite applies to LUMO and LUMO+1, which are thus localized on the functionalized 
part. While the numbers for IP and EA in [2] are taken from total energy calculations, therefore 
giving us good precision, we remark that the energy level spectrum shown here is by construction 
based on one-electron mean-field approximation and cannot give us a level-by-level accuracy, in 
particular for the unoccupied states. We see however that the overall picture is totally consistent 
with the type II alignment. 
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Conclusions 

In conclusion, we have studied the effects of organic covalent functionalization on sub-nanometer 
sized graphene nanofiakes, by varying the anchored groups, the number of functionalizing groups 
and the flake width. Our results indicate that edge functionalization can be designed to tune the 
electronic properties of graphene nanostructures and to control band alignment in graphene nano- 
junctions, while preserving stability and intrinsic electronic properties. This opens the possibility 
to realize stable all-graphene electronic and optoelectronic nanodevices based on type II offset.^ 
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